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1. Introduction 

The rare decay — H^jJ.^ is a process that is potentially sensitive to physics beyond the 
standard model (SM). In the SM, the decay can go only through penguin or box topologies at the 
loop level. Thus, a small branching fraction has been predicted, with the aid of lattice QCD, to be 
3.2(2) X 10"^ [|l|, |]. Recendy, LHCb [|] and CDF [Q] reported bounds on the branching fraction, 
to be followed by upcoming results from CMS. It is likely that a 5a measurement will be made, 
even at the SM branching ratio, in the near future. 

At LHCb, the extraction of the branching fraction relies on the normalization channels — )• 
J/yK^, B'^ K+7i^ and J/y^ [|], through the following relation 



BR(S0 ^ mV") = BR(B, ^ X)^^^, (1.1) 



where e and are the detector efficiencies and number of events. The fragmentation fractions, 
fq (q = u,d,s or A), denote the probability of a. b quark hadronizing into a Bq meson or a b- 
flavored (e.g., A^,) baryon. The fragmentation fraction ratio fs/fd is crucial in the extraction of 
BR(B^ — )■ n^n^). Currently, the uncertainty in fs/fd is the major source of uncertainty. Tra- 
ditionally, fs/fd was measured using the ratio of the corresponding semileptonic decays. Fleis- 
cher, Serra and Tuning proposed ||^ that the ratio can also be measured using the non-leptonic 
decays B^. — )• D^n and B*) — )• D^K , which has the advantages of a cleaner background, similar 
reconstruction of final states, etc. The approach is based on factorization of the nonleptonic am- 
plitudes into fji or fx and corresponding semileptonic form factors. The ratio fs/fd is related to 
BR(Bs — )■ Dn)/BR{B — DK) in a way similar to Eq. (1.1). With the efficiencies and event counts 
combined with the factorization approximation, we have 



fd \VudVfl T^o 



^ (1.2) 



where T is the lifetime and ^ ^ 1 with corrections of a few percent due to nonfactorizable ef- 
fects The semileptonic form factor ratio .y^f = [/d^'(M^)//o''^(M|-)]^ is currently the decisive 
contributor to the theoretical error. The value currently used at LHC^ is an estimate from QCD 
sum rules, = 1.24(8) ^. However, this theoretical input and the size of its error need to 
be validated by a nonperturbative method such as lattice QCD. This paper is devoted to such a 
calculation. 

The matrix elements of the B — )■ D semileptonic decay (and similarly for B^ — )• D,) can be 
written as 



+ (1.3, 

r 

However, for heavy quarks it is convenient to use the variables h±, defined by 

where v = p/M and the recoil variable is w = v • v' . We will use the form factors h± in our entire 



analysis and convert them in the end to /+,/o using Eqs. (1.3) and (1.4). 

In these proceedings, we report a preliminary result of the form factor ratio fjf^ (M^) /f^^^ (^1) 
by analyzing the semileptonic decays B^ — )- Dfl^v and B" — ^ D^l v on the lattice. We use an 
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identical subset of the MILC gauge configurations for both of the — )• Z), and B ^ D processes. 
To reduce the statistical errors effectively, we construct a set of ratios at small recoil, from which 
we extract the lattice form factors h±. The extrapolation to physical light quark masses and to the 
continuum is performed using root staggered chiral perturbation theory (rS;^PT). The results are 
extrapolated to maximum recoil by employing a model-independent parametrization. In Sec.^ we 
report our lattice result. 

2. Numerical details 
2.1 Data setup 

Our calculation uses four ensembles of the MILC's (2+l)-flavor gauge configurations two 
at each of the lattice spacings a f« 0. 12 fm and 0.09 fm. The ensembles as well as the parameters 
used are summarized in Table [T[ The strange and light sea quarks were simulated using the asqtad- 
improved staggered action [^. The action is also used in our strange and light valence quarks. The 
heavy quarks (charm and bottom) are simulated using the Sheikholeslami-Wohlert (SW) clover 



action with the Fermilab interpretation [10]. For the B — )• D decay, the spectator light quark is 
degenerate with the light sea quark (full QCD). While for the Bs — )• Dg decay, the strange quark is 
set close to its physical value. The charm and bottom quarks in our calculation are tuned to their 
physical values up to a tuning uncertainty. The corresponding bare hopping parameter K"^(c), as 
well as the coefficient for the clover term csw given explicitly in Table [l[ 

2.2 Lattice extraction 

In this work, we are interested only in the vector current operator. On the lattice we define 
V = y^'z^^Z^ mciy^^b, where Z'''' are normalization factors. The vector current in the continuum 

<,,Uq^q r\2 — 7bc 7cb /'-, 



is = pynV^, where Py^ = Zy^Z^p. jZ^nZ^a,. The factor can be calculated perturbatively and 



has been found to be very close to one [11, 12 1. We expect the pys to largely cancel in the ratio of 
the form factors. Hence this correction is negligible compared with other systematic errors and we 
take pv = 1 for this analysis. 

We employ the following three-point functions in our analysis, 

Cy^(0,f,r;p^) = £ (0|^z)(0,0)^,.r»F,(?,y)^+(0,x)|0) e*"^ (2.1) 
Cf;''^(0,r,r;pz,) = £ (0|^z)(0,0)^,/7''*I'.(^yX(0,x)|0) e'P"^ (2.2) 

Cf;;^(0,f,r;0) = £ (0|^B(0,0)^,//»P,(f,y)^+(0,x)|0). (2.3) 

The B meson is at rest. To obtain the dependence of the form factors at small recoil w, we simulate 
the final state D meson at a few small momenta, i.e., p = 27r(l,0,0)/L , 27r(l, 1,0)/L, 27r(l, 1, 1)/L 



a(fm) anil jams A^confs Csw 'Q- amj,{B ^ D) am^ (g,. D,, ) 

7^0. \2 0.020/0.050 2052 1.525 0.1259 0.0918 0.020 0.0349 

«0.12 0.007/0.050 2110 1.530 0.1254 0.0901 0.007 0.0349 

f«0.09 0.0124/0.031 1996 1.473 0.1277 0.0982 0.0124 0.0261 

«0.09 0.0062/0.031 1931 1.476 0.1276 0.0979 0.0062 0.0261 



Table 1: MILC ensembles of configurations used in this analysis. 
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and 27r(2,0,0)/L. The correlation functions for D — )• D and B ^ B serve as normalization. The 
D ^ D correlation function with a non-zero final state momentum is used to extract the recoil w to 
alleviate the need of renormalizing the four velocity. 

From these correlation functions we construct three different ratios, fits to which include con- 
tributions from the lowest-order excited states. Explicitly, 



c 



3pt 



CZ'''{0,t,T;p) 



d' (1 + %2 + ^20 ^'-^^(P^O, 



-^DVB 



cg;'^(o,^^;p) 

Cg;^^(0,f,r;0) 



(2.4) 
(2.5) 



Zo(0) / £o(p) ^(g„(p^_j;„(o^,, 
Zo(p) Y £o(0) 

' (1 + ^02 e-'^'^^^-'^+^zo 



(2.6) 



The factor in the square brackets of Eq. (2.6) cancels the time dependence of the ratio, stemming 
from the fact that the numerator and denominator have final state D mesons with different momenta. 
AS and Am denote the lowest splittings and 5 is a parameter that accounts for the imprecise ^(p) — 
£■(0) in the bracket in Eq. (p^). In the fits the lowest-lying energy splittings AE, Am are treated 
as fit parameters. The splittings can be extracted from the two-point functions. So, we employ a 
multi-channel fitting procedure, combining the two-point functions and the ratios of the three-point 
functions. We find that such a treatment results in more robust fits and more precise splittings. 
From di,bi and a, we can easily recover the form factors h± at small recoil w, 

1+d d 



1-d d' 

h+{w) = /j+(l)(a//Z7i-a-d), 

h-{w) = h+{\){ai/bi- Qi/di). 



(2.7) 

(2.8) 
(2.9) 



3. Results 

The extrapolation of our lattice results to the physical quark masses and the continuum is 



guided by rS;tPT [13, 14]. However, in the case of h+ the light quark mass dependence is accom- 
panied by a small recoil w dependence. Such a dependence was included in the continuum chiral 
perturbation theory in Ref. [15] and was extended to the NLO rS;^PT in Ref. [^]. For the 
NLO correction is simply a constant which is inversely proportional to the charm quark mass. We 
follow the same setup, adding NNLO analytic terms and including dependence. The remaining 
recoil dependence of the form factors is fitted to a simple quadratic expansion at zero recoil. 

The results of the chiral/continuum extrapolation are shown in Fig. |l|. The form factor /j+ for 
both of the B — )• D and Bg — )• D., decays shows a small dependence on the light quark masses and 
lattice spacings. The extrapolated physical values are very close to the lattice data points. This 
suggests that /j+ is insensitive to the light degrees of freedom. However, sizable light quark mass 
and lattice spacing dependence appears in the case of , as indicated by the variation due to the 
sea quark masses and the differences between and /j^'^^' (spectator mass). Note that the 

difference between and h^^^' is minor Considering the subleading role that /j_ plays in 

contributing to f+,fo, we expect the U-spin symmetry breaking effect to be smaller than what was 
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expected in ||^ 0]. Such an observation is bolstered by the recent lattice calculations on /+,/o of 
the D(,) n{K) decays 

With the physical values of h±, we can easily calculate /o,/+ using the physical masses of 
the D and B mesons. However, to evaluate the form factors at a small momentum transfer (q^ = 
M^,M^), we need to extrapolate the results near maximum recoil. We use the model-independent 



z-parametrization [18] with the consti^aint /o(0) = /+(0). We take four synthetic points in the 
recoil range where we have lattice data points. We take the values of /+,/o by evaluating our 
chiral/continuum extrapolation result at these four recoil points and perform the z-parametrization. 
The result is shown in Figure ^. We study the effect of a pole at a vector meson in the Blascke 
factor of the z-expansion of /+. We find that the shapes of the form factors are only weakly affected 
by the inclusion of such a pole. 

By expanding the form factors at the respective momentum transfers, we finally arrive at 

fl^\Ml)/i'\Ml) = 1.05 1(40) (22). (3.1) 



The first error is from statistics. The second error is the systematical error due only to the un- 
certainty on gDD'n and to the variation of fits in the z-parametrization. We are in the process of 
building a full systematic error budget. 
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B->D z-expansion: x"/dof=0.48/l, CL=0.49 B,->D, z expansion: x^/dof=1.3/l, CL=0.25 




Figure 2: The z-expansion of form factors /o,/+. The points that we include in the z-expansion fits are 
shown expHcitly. The dashed curves indicate the result of chiral/continuum extrapolation. 

4. Conclusions 

In summary, we present a (2-1-1 )-flavor lattice QCD calculation of the form factor ratio 
/q*^(M^)//q''^(M|^), which is a major theoretical input for the extraction of the fragmentation 
fraction ratio fj fd- The essential part of our calculation is to extract the B ^ D and Bs — )• 
semileptonic form factors at non-zero recoil. We reduce the systematic uncertainty by fitting the 
lowest-order excited states, and we employ a simultaneous multi-channel fit procedure to address 
correlations and reduce the statistical uncertainty. Our chiral/continuum results show that the cor- 
rections to the finite lattice spacings and finite light quark masses are small. Our preliminary result 
is /o^^(M^)//d'^' {^k) ~ 1-05 1(40) (22), with a partial systematic error budget. As a consequence, 
we obtain jVp = 1.105 (80) (44) which implies a smaller U-spin breaking effect than that suggested 
in [^, c/ff = 1.24(8). A more comprehensive analysis with a detailed error budget is still in 
progress and will be reported in a forthcoming paper. 
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